Evidence of a glass transition in a 10-state nonmean-field Potts glass 
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Potts glasses are prototype models that have been used to understand the structural glass tran- 
sition. However, in finite space dimensions a glass transition remains to be detected in the 10-state 
Potts glass. Using a one-dimensional model with long-range power-law interactions we present ev- 
idence that a glass transition below the upper critical dimension can exist for short-range systems 
at low enough temperatures. Gaining insights into the structural glass transition for short-range 
systems using spin models is thus potentially possible, yet difficult. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 



Due to their fascinating properties such as aging, mem- 
ory effects and ergodicity-breaking transitions, as well 
as industrial applications, structural glasses, supercooled 
liquids and polymers have received considerable atten- 
tion recently. In particular, when the temperature is de- 
creased, they undergo a dynamic transition [H-IH below 
which the particle-density correlation length does not de- 
cay to zero in the long-time limit and the evolution be- 
comes noncrgodic. However, this transition is not asso- 
ciated with any thermodynamic singularity. Hence the 
system "freezes" in a portion of phase space. There is 
a second transition at a lower temperature [1, Q which 
can be associated with a thermodynamic singularity and 
which can be related to a possible ideal glass transition. 
Despite ongoing efforts, the structural glass transition 
remains to be fully understood. 

The p-state Potts glass @4LH is one of the most ver- 
satile models in statistical physics: For p = 2 states it 
reduces to the well-known Edwards- Anderson Ising spin 
glass [lH, a workhorse in the study of disordered mag- 
netic systems. For p = 3 it can be used to model orienta- 
tional glasses [l3[ , while for p = 4 the Potts glass can be 
used to model quadrupolar glasses. For large p > 4 and 
no disorder the model shows a first-order transition. In 
particular, infinite-range Potts glasses with p > 4 exhibit 
a transition from ergodic to nonergodic behavior [r34TT|]. 
as well as an additional static transition at a lower tem- 
perature. In fact, the equations describing the system's 
dynamics near the transition are mathematically related 
[14H16l ] to the equations of mode-coupling theory, which 
describe the behavior found in structural glasses and su- 
percooled liquids. Therefore, studying the Potts glass 
with large p could provide, in principle, some insights into 
the mechanisms governing the structural glass transition. 
However, this beneficial relationship seems to only work 
when the model is infinite ranged 17]. The existence of 
a transition in finite-dimensional systems remains to be 
proven [l8|, [l9[ . Not only arc hypcrcubic lattices with 
large space dimension hard to study numerically, recent 



work [20j,|21| suggests that if there is a transition for large 
p it would occur at very low temperatures. 

In this work we simulate the 10-states Potts glass on 
a one-dimensional ring topology with power-law interac- 
tions. This allows us to effectively tune the range of the 
interactions and therefore the (effective) space dimen- 
sion for large linear system sizes. Our results suggest 
that 10-state Potts glasses should have a very low finite- 
temperature transition for finite space dimensions. 

The paper is structured as follows. In Sec.|H]we intro- 
duce the model and observables. Furthermore, we outline 
the details of the numerical simulations. Section |HT] sum- 
marizes our findings, followed by concluding remarks. 



II. MODEL AND OBSERVABLES 

We study a one-dimensional Potts glass with long- 
range power-law interactions (22l . |23| and Hamiltonian 
U = - J2 itJ -Jijfiqiw, where q l € {1, ...,10} are 10- 
state Potts spins on a ring of length L to enforce pe- 
riodic boundary conditions and 8 x ^ y = 1 if x = y and 
zero otherwise. The sum is over all spins and the in- 
teractions Jij are given by Jij = Eij/rfj, where are 
Normal distributed with mean Jo and standard devia- 
tion unity. Tij = (L/tt) sin[(7r|i — j\)/L] represents the 
geometric distance between the spins on the ring. For 
the simulations we express the Potts glass Hamiltonian 
using the simplex representation where the 10 states of 
the Potts spins are mapped to the corners of a hypcrte- 
trahedron in nine space dimensions. The state of each 
spin is therefore represented by a nine-dimensional unit 
vector Si taking one of the 10 possible values satisfy- 
ing the condition ■ S v = [p/(p - 1)](<W - !) with 
{//, v} G {1,2,..., 10}. In this representation the Potts 
glass Hamiltonian is given by % — — ■ Jij Si ■ Sj with 

Jij = Jij(p — l)/p. In the limit when a — > 0, when the 
system is infinite ranged (Sherrington-Kirkpatrick limit), 
we obtain T c {a = 0) = l/(p - 1). 

The merit of the long-range one-dimensional model 
lies in emulating a short-range topology of varying di- 
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mcnsionality, depending on the power-law exponent: For 
a < 2/3 the model is in the mean- field long-range 
10-state Potts universality class and, in particular for 
a < 1/2 in the infinite-range universality class. However, 
for 2/3 < a < 1 the model is in a nonmean-field uni- 
versality class with a finite transition temperature T c . It 
can be shown [22[ that a = 2/3 corresponds exactly to six 
space dimensions for a hypcrcubic lattice. Therefore, a 
values between 1/2 and 2/3 allow us to effectively study 
}22l l23j a short-range hypercubic Potts glass above the 
upper critical dimension d u = 6, whereas when a > 2/3 
we effectively study a model with a space dimension below 
six dimensions. Thus, by studying the one-dimensional 
model we can infer if a transition should be present for 
the corresponding short-range hypercubic Potts glass. 

The presence of a transition is probe d by studying the 
two-point finite-size correlation length [24[ • We measure 
the wave- vector-dependent spin-glass susceptibility (25j 



XsG (k)=iV^[(|^(k)| 2 )] £ 



where (■ • • ) denotes a thermal average, 
over the disorder and 



(1) 



an average 



(k) = ^£sf (Q) ^e^ 



(2) 



is the spin-glass order parameter computed over two 
replicas (a) and (/3) with the same disorder. The two- 
point finite-size correlation length is then given by 
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2sin(/c min /2) 



Xsg(O) 

XSG(k m in 



1/(2,7-1) 



(3) 



where k m j n = 2-k / L is the smallest nonzero wave vector. 
According to finite-size scaling }25j 

a/L v/3 = X[L 1 / 3 (T-T c )] (1/2 < a < 2/3) , (4a) 
fi/i - X[L 1/lJ {T - T c )] (2/3 < a), (4b) 



TABLE I: Parameters of the simulations for different ex- 
ponents a. N sa is the number of samples, iV sw is the total 
number of Monte Carlo sweeps, Tmin is the lowest tempera- 
ture simulated, and Nt is the number of temperatures used 
in the parallel tempering method for each system size L. 
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[Eq. (P) with k = 0]. 
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In general, the exponent rj has to be known a priori to 
precisely determine the location of T c . However, for the 
one-dimensional model 2 — ?/ = 2a — 1 for a > 2/3 ex- 
actly and so XSG/L 2 ~ n can be treated as a dimensionlcss 
quantity similar to the two-point correlation length. 

To prevent ferromagnetic order 0, 0] we set the mean 
of the random interactions to Jo = — 1 fl8j in our simu- 
lations. This suppresses the ferromagnetic susceptibility 

Xrn = iV£ M [(l™f }]av [m" = (l/JV) 5?]. We disCUSS 

the case where Jo = — 1 in more detail below. 

The simulations are done using the parallel temper- 
ing Monte Carlo technique prjj ; simulation parameters 
are shown in Table [T] Equilibration is tested by using 
an exact relationship between the energy and four-spin 
correlators (link overlap) [27j when the bond disorder is 
Gaussian, suitably generalized to Potts spins [l9| on a 
one-dimensional topology [28| . 



where v is the critical exponent for the correlation length 
and T c the critical temperature. For a < 2/3, v = 

1/(2(7-1). 

In practice, there are corrections to scaling to Eqs. ((4]) 
and so data for different system sizes do not cross ex- 
actly at one point as implied by the finite-size scaling 
expressions. The crossings between pairs of system sizes 
L and 2L shift with temperature and tend to a con- 
stant for L -» 00. In general, T* = T r °° + b/L 8 with 
9 = 1/v + lo. Here we find empirically that l/f + ufsl. 
We fit T*(L,2L) with high probability to a linear func- 
tion in 1/J. The intercept with the vertical axis after the 
fit determines a lower bound for the transition tempera- 
ture. Error bars are determined via a bootstrap analysis. 

To obtain a better understanding of the corrections 
to scaling we also measure the spin-glass susceptibility 



III. RESULTS 

Our results are summarized in Fig. [1] The main panels 
in the left column show data for the finite-size correla- 
tion length as a function of temperature for (a) a = 0.60, 
(c) 0.75, and (e) 0.85. The insets show the correspond- 
ing data for the scaled dimensionlcss susceptibility. In 
all cases data for different system sizes cross, indicat- 
ing the presence of a transition. To better quantify the 
thermodynamic behavior, we show in the right column 
the scaling of the crossing between successive system size 
pairs T*(L,2L) as a function of 1/L. The data can be 
well fit by a linear function; the intercept with the vertical 
axis corresponding to the thermodynamic limit. For all a 
studied we find finite values for the thermodynamic glass 
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FIG. 1: (Color online) Panels (a), (c) and (e) show the correlation length £l/L (inset: susceptibility xsg/L ~ v ) as a function 
of temperature T for different system sizes L. Panel (a) shows data for a — 0.60 (mean-field regime) where a transition is 
expected pffll ] [note that here v = l/(2<r — 1)]. Panels (c) and (e) show data for a = 0.75 and a = 0.85, respectively, which 
correspond to a space dimension below the upper critical dimension. A transition for low yet finite temperature is clearly visible. 
Panels (b), (d) and (f) show the crossing temperatures T*(L,2L) of successive pairs of system sizes for different exponents 
a [(b) 0.60; (d) 0.75; (f) 0.85]. The crossings for both £,l/L and Xsg/L 2-71 are well approximated by a linear behavior in 
1/L. Despite small deviations between the estimates for both quantities, for all a values studied T c (a) > 0. In particular, we 
estimate T c (0.60) = 0.060(4), T c (0.75) = 0.040(3) and T c (0.85) = 0.025(3). Note that the data for a = 0.60 show a deviation 
from the linear behavior for the largest system sizes studies. However, both data sets agree and therefore suggest that the 
thermodynamic limit might have been reached. 



transition. These findings for the long-range model with 
power-law interactions imply that the 10-state mean-field 
Potts glass, for d u < d < oo space dimensions, has a sta- 
ble glass phase at finite temperatures. In addition, our 
data for cr > 2/3 indicate that short-range Potts glasses 
with a space dimension below the upper critical dimen- 
sion should also have a finite transition temperature, al- 
beit at very low T [30|. 

Recently, Alvarez Banos et al. [2l[ performed a thor- 
ough study of a three-dimensional Potts glass with p < 6, 
bimodal disorder and Jq = 0. Their main result is that T c 



decreases with an increasing number of states p and sug- 
gests that for 10 states T c should be strongly suppressed, 
in agreement with our results. In addition, Alvarez Bahos 
et al. [2l| claim that (1) only weak ferromagnetic order 
is visible when Jo = 0, (2) that the complexity of the 
simulations is much higher when Jo = 0, (3) that setting 
Jo = —1 could impact the presence of the glass transi- 
tion, and (4) that the transition could be first order. 

We have examined these claims using the one- 
dimensional model with Gaussian disorder and find that 
(1) ferromagnetic order grows considerably when Jq = 
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FIG. 2: (Color online) Ferromagnetic susceptibility \ m as a 
function of temperature T for different system sizes. Ferro- 
magnetic order is strongly suppressed for Jo = —1 in compar- 
ison to the Jo = case. 

at low enough temperatures (see Fig. [2]) and (2) the com- 
plexity of the simulations is not affected by shifting the 
mean of the interactions. With respect to point (3), we 
do find, however, that the transition temperatures are re- 
duced by approximately a factor of 2 - 3 when Jo = — 1 
in comparison to the simulations where Jq = 0. Shift- 
ing the mean of the interactions therefore only quantita- 
tively impacts the transition temperature. Finally, (4), 



for the system sizes studied, the distribution functions of 
the energy show no double-peak structure that would be 
indicative of a first-order transition. 



IV. CONCLUSIONS 

Using a one-dimensional 10-state Potts glass with 
power law interactions, we present evidence suggesting 
that short-range finite-dimensional 10-state Potts glasses 
should exhibit a finite-temperature transition for low 
enough temperatures and large enough system sizes. Al- 
though corrections to scaling are large, we estimate that 
for all a values studied T c (cr) > 0. In particular, we 
conservatively estimate T c (0.60) = 0.060(4), T c (0.75) = 
0.040(3), and T c (0.85) = 0.025(3). Larger system sizes 
might show a different behavior, however, the presented 
state-of-the-art simulations show strong evidence that 
short-range 10-state Potts glasses in high enough space 
dimensions should order. 
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